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We consider a chain of free electrons with periodically switched dimerization and study 
the entanglement entropy of a segment with the remainder of the system. We show that it 
evolves in a stepwise manner towards a value proportional to the length of the segment and 
displays in general slow oscillations. For particular quench periods and full dimerization an 
explicit solution is given. Relations to equilibrium lattice models are pointed out. 

I. INTRODUCTION 

Entanglement is a basic feature of non-trivial quantum states and has been studied extensively 
for the ground states of standard many-body systems However, the question how it evolves in 
time-dependent situations is equally interesting. The simplest case is a quench where the Hamil- 
tonian is changed instantaneously from Hq to H\ and one follows the subsequent evolution of the 
initial state, usually taken to be the ground state of Hq. For global quenches, where a parameter is 
changed everywhere in the same way, it was found that in homogeneous chains the entanglement 
entropy S between a segment of length L and the remainder first increases linearly with time 
and then approaches an asymptotic value proportional to L . Thus the time evolution 

leads to a much larger entanglement than found in equilibrium where even for critical chains S 
is only proportional to InL. This result can be interpreted in a simple picture where pairs of 
quasiparticles establish the entanglement between the two parts of the system j^j]. Alternatively, 
the time-dependent state may be viewed as the ground state of an evolving Hamiltonian which is 
strongly non-local and gives long-ranged correlations. 

In this paper we will study a more general situation and consider not a single quench, but a 
periodic sequence of changes Hq *-* Hi and its effect on the entanglement entropy. Such periodic 
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changes have been considered in some papers on Ising chains kicked by a transverse 
or a tilted field B fl Q. We will investigate a system of free electrons hopping on a chain 
with alternating bonds. The corresponding dimerization parameter is then switched between the 
values ±5, which simply interchanges weak and strong bonds. Although there is a mapping to the 
transverse Ising chain, the hopping model has a more direct physical significance since it could be 
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realized with optical lattices. There have already been experiments with dilute gases where periodic 
potentials were applied in the form of kicks . Other experiments [ijl, [3] have motivated studies 
of chains where the parameters in the Hamiltonian were assumed to vary harmonically. In this 
case, however, the focus was either on the energy absorption Q, [lfj or on work fluctuations 



1 

The periodic quench has some important general features. In the evolution, one can distinguish 
what happens after full quench periods and what happens at the times in between. This situation is 
analogous to the case of a particle moving in a periodic potential and has been discussed extensively 
for simple kicked quantum systems 



181 ] . For a full quench period, the time evolution operator is 
the product of two exponentials and analogous to the transfer matrix in a classical two-dimensional 
system with a layered structure. One can write it as a single exponential with a certain effective 
Hamiltonian. The formulae are very similar to those one finds for the Ising model on a square 
lattice. They also show that the effective Hamiltonian is critical if, as we assume, Hq and H\ act 
for the same length of time. Compared to a simple critical quench, however, one has the additional 
evolution during the intermediate times. The unitarity of the discrete evolution from period to 
period, on the other hand, leads to effects not present in models with Euclidian time. 

We describe the model and the method of computing S via correlation functions in Section 
2. Then we treat in Section 3 the fully dimerized case 5 = 1 which is particularly simple. The 
time evolution can then be represented graphically and for special values of the quench period 
the particles even move with strictly uniform velocity. For S(t) one finds a step structure and an 
asymptotic value proportional to L. We discuss the dependence on the quench period and show the 
behaviour for arbitrary times. Section 4 deals with arbitrary dimerization. Here the quasiparticles 
in the effective Hamiltonian can have a more complicated dispersion in which case additional slow 
oscillations appear in the entanglement entropy. In Section 5 we sum up our findings and in the 
Appendix we give the formulae for calculating the correlation matrix in the general case. 

II. MODEL AND METHOD 

We will study a half-filled ring of free electrons with 2N sites and alternating hopping as shown 
in Fig. [TJ 

The Hamiltonian is 
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H = -J ( C 2™ C2m+1 + C 2m+l C 2m+2 + h.C. J (1) 

m=l ^ ' 
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FIG. 1: Geometry of the dimerized hopping model. 

and, up to boundary effects, also describes a dimerized XX spin model. In the following we set 
J = 1. The diagonal form of H is obtained by working with cells of two sites, defining C2 m = a m , 
C2m+i = b m , and performing a Fourier transformation 
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Then 



ff = -j> g (aja, (3) 
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where 



cos 2 | + 5 2 sin 2 | (4) 



and the Fermi operators are related via 



with 



_L e -w 2 K + Az ) ; ^ = -^6^(0,-^) (5) 



,v,/9 cos o — i^sin % , v 

e^ 9 = e 4g/2 ^ 2- . 6 

In the ground state, all a-levels are occupied and all /3-levels are empty. The single-particle 
spectrum has a gap at q = ±7r for all 5^0. For vanishing dimerization 5, (JH) gives u q = cos(g/2) 
which is the usual result for homogeneous hopping, but folded into the smaller Brillouin zone. For 
5 = 1 the dispersion becomes flat, uj„ = 1. A change 5 — > —5 does not affect co q which simplifies 
the situation in the quench. 

At time t = 0, the system is taken to be in the ground state l^o) of -ffo = H($). The periodic 
quench consists of switches between Hq and H\ = H{—5) at times nr where n = 0,1,2, ... . This 
interchanges weak and strong bonds and forces the system to adapt. Formally, the wave function 
then evolves according to 

|<&(t)> = U(t)\$ ) (7) 



with a unitary operator U(t) which for t = n ■ 2t, i.e. after n full periods, is given by 

U(2nr) = U n (8) 

where 

U = U U 1 = e - iHoT e- iHlT (9) 

describes the evolution over one period. 

We are interested in the entanglement of a segment of L sites with the remainder of the sys- 



tem. The corresponding entang" 



segment which has the form 
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ement entropy S follows from the reduced density matrix p of the 
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Here Z is a normalization constant ensuring Tr p = 1 and the fermionic operators fk follow from 
the c n by an orthogonal transformation. Then S = — Tr {p\o.p) is determined by the single-particle 
eigenvalues according to 

= In Cfc(t) " E^ 1 " Cfc(*))Ml " Cfc(*)) (ii) 

where Cfc(i) = l/( ex P( £ fc(*)) + !)• The Cfc(^) are the eigenvalues of the correlation matrix 

C lm (t) = (<5>(t)\4c m \<S>(t)) = (<S> \4(t)c m (t)\<S> ) (12) 

restricted to the sites of the subsystem. To obtain C(t) one relates it to C(0) by expressing the 
Heisenberg operators c m (t) = W (t)c m U(t) as linear combinations of the Cj. For this, one has to 
diagonalize the operator U in ([U]) which is the product of two exponentials. This is a well-known 
problem which first appeared (without the factors of i) in the solution of the two-dimensional 



Ising model via transfer matrices 2(| • A special case is treated in Section IIIII while the general 
expressions for the elements of C(t) are given in the Appendix. Evaluating them numerically and 
diagonalizing the matrix one then finds the entropy. 

One should note that there is a close relation between XY and transverse Ising (TI) chains 



21 



22l | which was recently extended to their entanglement properties This means that 

the Hamiltonian H in (pQ) can be obtained by putting two TI models on alternating sites of a 
chain, performing a dual transformation and choosing the coupling constants \ a and the fields h a 
as h\ = A2 = (1 + <5)/2, h% = Ai = (1 — 5)/2. Thus one TI model is in the ordered and the other 
in the disordered phase. Changing the sign of 5 then switches each model to the opposite phase. 



In the special case 5 = ±1 one has either only a field or only a nearest neighbour coupling. This 
corresponds to an Ising chain kicked periodically by a transverse field {(], 7, 8], because during 
the 5-pulses describing the kicks one can neglect the coupling terms. One could work in this 
TI formulation, but the dimerized chain provides a nicer physical picture as will be seen in the 
following section where we actually treat the case 5 = 1 before dealing with the general situation. 



III. FULLY DIMERIZED CASE 



In this section, we treat a periodic quench between 6 = ±1 which is particularly transparent 
and instructive. This situation has also been considered in 25] for bosons. 



A. Time evolution operator 



For 5 = 1, the system consists of independent pairs of sites. For each pair, the Hamiltonian has 
the form 



h 



and is diagonalized by setting a = (a + (5)/^/2, b = (a — (3)/^/2. This corresponds to 



if q = in 



and gives 
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b 



(14) 



The time dependence of a and b is then 
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cos t i sin t 
i sin t cos t 

and the full evolution matrix has block form with N such v's along the diagonal. The same holds 
for the other half-period but then the pairs and thus the v's are shifted by one site. The complete 
time evolution can therefore be represented as in Fig. [51 where each shaded square corresponds 
to a matrix v(t). Such a chequerboard structure also appears in the time evolution of lattice- 
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gas models with parallel dynamics 
functions of vertex models oriented diagonal 
with a Trotter decomposition of exp(—f3H) 



281 ]. In equilibrium problems it is found for partition 



y 2j| or of one-dimensional quantum systems treated 



Using translational invariance, the (2iV x 2iV) matrix V = VoVi for one period is found to 
have the eigenvalues exp(zkij q ) where 

(16) 



cos 7 g = cos 2 r — sin 2 r cos q . 
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FIG. 2: Time-evolution pattern for full dimerization. The shaded squares indicate the coupling of the sites 
and of the corresponding Fermi operators in each half-period. 



Therefore, writing the operator U in ([9]) as a single exponent 



U = e 



the diagonal form of the effective Hamiltonian is 



-iH2j 



(17) 



(18) 



where v q = j q /2r and the Fermi operators £ 9 , r/g follow from the eigenvectors of V. 

The single-particle energies of H are shown in Fig. [3] for several values of r. The main feature 




FIG. 3: Dispersion of the single-particle excitations in H for five values of the half-period r. 



is that the spectrum is always gapless at q = ±7r. Thus the evolution over one period is given by a 
critical Hamiltonian. This is a consequence of the symmetry of the two half-periods. If they have 
different lengths tq, tl, Eq. (fT6|) is changed to 



cos 7 g = cos tq cos ti — sin tq sin t\ cos q 



(19) 



and 7^ = |to — T\ \ . One should note that the relation (fT9|) is the trigonometric analogue of the 
hyperbolic formula 

cosh j q = cosh 2K * cosh 2K2 — sinh sinh 2K2 cos g (20) 

found for the row transfer matrix of the two-dimensional Ising model on a square lattice with 
couplings K\ and K2 The quantity K\ in (I20p is the dual coupling of K\, tanhiC* = 

exp(— 2K\). This illustrates the close connection between the two problems mentioned earlier. 



B. Hamiltonian limit 



For rapid switching, r<l, one finds from (|16|) j q = 2tcos |, i.e. v q = cos |. According to (JH) 
these are the single-particle energies u; 9 of H in the absence of dimerization. This is not surprising 
since for t<1 one can write directly 

U = e - iH ° T e~ iHlT ps e-^o+^iK (21) 

which gives -ff = (-ffo + Hi) /2 = H(5 = 0). Therefore the evolution for times t = n ■ 2r is the same 
as for a single quench to a homogeneous hopping model. This situation has already been studied in 



311 ]. Appendix B, where the correlation matrix C(t) was calculated and the entanglement entropy 



was found to approach the asymptotic value 

5(oo) = L(21n2 - 1) (22) 

for large L. 



C. The special case r = tt/2 

For this value of r, v q = 1 ± q/ir is strictly linear and describes massless particles with velocity 
v = 2/tt in units of the lattice spacing. This can also be seen directly in real space. According to 
(fT5j) . one has 

aC-) = ib , 6(|)=ia (23) 

such that particles in a cell interchange positions. This continues in subsequent half-periods and 
leads to straight paths as sketched in Fig. HI 

A correlation which exists initially between the first neighbours A and B is thereby extended 
to third neighbours in the first step, to fifth neighbours in the second one and so forth. This 
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FIG. 4: Motion of the particles for r = 7r/2. 

permits to set up the correlation matrix C(t). If A and B form a cell, then in the ground state all 
correlation functions are equal 

(ata) = (atft) = (&t ) = (fttft) = I (24) 

and the same holds then for correlated sites at later times. For the entanglement entropy only 
those correlations are relevant which connect sites inside and outside the subsystem. They lead to 
an eigenvalue £ = 1/2 which gives a contribution of In 2 to S. It is not difficult to see that in one 
step the number of in-out correlations does not change while in the other one it increases by 4. 
Thus in a full period 

AS = 4 In 2 (25) 

and the increase continues until S reaches the highest possible value Lin 2. This build-up of the 
entanglement is an exact lattice realization of the light-cone picture introduced by Calabrese and 



Cardy [2] in the context of single quenches and analyzed further in 32. l33|. 

The description becomes complete if one also looks at intermediate times. By studying small 
values of L one finds that during the half-period where S changes, the correlation matrix has 
non-trivial eigenvalues £i = sin 2 (t/2) and £2 = cos 2 (i/2). With (llip these give a contribution 



AS(t) 



sin 2 — In sin 2 — h cos 2 — In cos 2 — 
2 2 2 2 



(26) 



which describes the increase by 4 In 2 between t = and t = tt/2. Therefore S(t) has a step-like 
structure where plateaus are connected by the function AS(t). The result of a numerical calculation 
for L = 20 is shown in Fig. [5l Due to the alternating structure of the chain, there are two possible 
choices for the subsystem if L is even. It may contain initially L/2 complete cells, in which case 
it is uncorrelated with the surrounding and S(0) = 0. Or it contains only half a cell at each end, 
which contributes In 2 to the entanglement giving altogether S(0) = 2 In 2. This leads to the two 
curves in the figure. Note that for the second choice the last plateau is absent and the final step 



has only height 2 In 2. Its form is also slightly modified. The inset shows a comparison of AS(t) 
with the analytical formula ([26]). 
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FIG. 5: Entanglement entropy S(t) for r = 7r/2, L = 20 and the two choices of the subsystem, see text. 
Inset : Enlarged view of the ascending part AS(t). The line is the analytical result. 

These considerations can also be extended immediately to times r which are multiples of 7r/2. 
Then the particles change their position in the cells several times during one half-period. Conse- 
quently, the plateaus are longer and S(t) follows the function (|26|) for a longer interval. Thus if 
r = (2fc + l)7r/2 the plateaus are connected by an oscillating function. Such a case is also discussed 
below for more general times. On the other hand, if r is an even multiple of vr/2, the oscillations 
always lead back to the initial value of S and there is no systematic ascent. 



D. General times r 



For general values of the quench period, no analytical treatment is possible and the entropy has 
to be obtained numerically as described in Section 2. Figure [6] shows results for small values of r 
on the left and for larger ones on the right. 

The overall behaviour is an initial, on average linear increase followed by a sharp bend and a 
final approach to an asymptotic value. The upper bound Lin 2 is only reached for r = 7r/2. The 
bending points can be related to the velocity of the fastest one-particle excitations in the effective 
Hamiltonian. This is given by 



^max — 2 



dq 



(27) 
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FIG. 6: Entanglement entropy S(t) for L = 40 and different values of r. Left: Small r. Right: Larger r. 
The arrows indicate the times L/2i> max . 



where the factor of 2 is needed for the distance to be measured in sites instead of cells. Then 
the time required for these excitation to travel from the middle of the block to the surrounding 
is given by L/2u max . The corresponding values are indicated by arrows in Fig. [6] and locate the 
cross-over points very well. In the Hamiltonian limit r — > 0, the result for the single quench to the 
homogeneous chain is recovered. 

Let us now discuss the fine structure of the curves. First of all, the increase takes place in steps 
with well defined plateaus. These plateaus still occur because in one of the two half-periods the 
time evolution in the subsystem is decoupled from that in the surrounding. Thus the entanglement 
cannot change. For small r, the plateaus are connected by simple monotonous curves, but for 
larger periods, an overshooting occurs which turns into an oscillation for even larger r. A particle 
then performs more than one cycle in a cell before the next switching occurs. This case is shown 
in more detail in Fig. [7|for r = 1.7ir. The first oscillation is given by the function (j26[) up to the 
time t = t when the plateau starts. The next oscillations look similar but their amplitudes are 
not the same and even after a rescaling there remain small differences. In fact, their similarity is 
even somewhat surprising. On the right of Fig. [7| the non-trivial eigenvalues which determine 
the variation of S(t) are plotted. One sees that only one such eigenvalue (and its partner 1 — 
exists in the beginning but that more and more appear in later periods and contribute to S. Thus 
one expects differences in the functional form. 

Looking at the step heights in Fig. [7J one sees that they are not all equal. This is a general 
feature and illustrated in Fig. [SJ The main graph shows the step height AS n of the n-th step for 
the initial increase of S(t) with the average value subtracted. It is an oscillatory decaying function 




4 
tlx 



2 4 
t/x 



FIG. 7: Left: Entanglement entropy for S(t) for L = 40 and r = 1.77T. Right: The eigenvalues Cfe(X) 
producing S(t). 



independent of L and can be very well described by the expression 

sin(4rn — ip) 



AS n = AS + A 
with an amplitude A and a phase <p. 
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FIG. 8: Variation of the step height in S(t) for L = 200 and r = 0.17T. The average value ASq is subtracted. 
Inset: A5o as a function of r. 

This result is easy to understand if one takes into account that the integrals (I39p for the elements 
of the correlation matrix contain the number of the period in the form of factors sin(27 g n) and 
cos(27 g n) and that the main contribution comes from the region near q = where the dispersion 
curve has its maximum. This gives oscillations in n with frequency 270 = 4r which persist in S(t). 
Moreover, the density of states in the integration leads to a dependence l/n 3//2 , which was also 
found in a previous treatment of the Hamiltonian limit 1 311 and in a calculation of a related TI 

n 

correlation function [g|. The oscillations go away if either r = 0orr = 7r/2. Between these two 



limits, the average step height increases monotonously from zero to 4 In 2 as shown in the inset of 
the Figure. Thus the case r = n/2 studied in the last subsection has not only regular steps but 
also the highest ones and thus leads to the fastest increase of S. Since the dispersion relation ^ q 
is invariant under r — > tt — r and under r — ► r + kir, one can use the above results to obtain the 
step heights for arbitrary periods. 

In the region beyond the cross-over point cit time L/^v^^ the step height decays rapidly towards 
zero with the same n~ 3//2 power law for large n as before and one can still observe the same 
oscillations. 

IV. PARTIALLY DIMERIZED CASE 

We now turn to the case of general values 5 < 1, where the system has alternating weak and 
strong hopping. Results for a short quench time, r = QAir, are shown in Fig. [9l The curves 
resemble those in Fig. \5[ but there is one basic difference. Since the cells are always coupled, one 
does not have plateaus any more. Instead, S(t) increases slowly during those half-periods where the 
plateaus occur for 5 = 1 and faster for the others. Together this still gives a kind of step structure, 
but the steps are washed out. As 5 becomes smaller, the differences between the two half-periods 
also become smaller and the curve smoothens more and more. At the same time the average slope 
decreases. Moreover, the initial value S(0) is non-zero and becomes larger, approaching the value 
S = 1/3 In L + k of the homogeneous chain. The oscillations with the quench period persist also 
beyond the bending point and decrease there with 5 in the same way. 
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FIG. 9: Entanglement entropy for L = 40 and r = OAw for five values of the dimerization 5. 



This might suggest that the behaviour of S becomes less interesting for partial dimerization. 
However, this is not true. In Fig. [TUlwe show results for a larger quench time, r = 0.77T and several 
values of 5 below 0.5. In addition to small and fast oscillations one now observes very slow ones 
with rather large amplitudes and a frequency which decreases with 5. As in the previous section, 
the origin of these oscillations can be found in the dispersion curve for j q . In Fig. [11] we have 
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FIG. 10: Entanglement entropy for L — 40 and r = 0.7tt for several values of the dimerization 5. The dotted 
lines on the right indicate the asymptotic values. 

plotted 7 g given by Eqs. (p?2"]) and in the Appendix as a function of S. One sees that while 
there is a single maximum at q = for large 5, a double-peak structure develops below 6 = 0.7 
and the maxima approach the zone boundary at tt for 5^0. The reason is that for small 5 the 
factor cos((/?i — ipo) is almost constant and -f q is basically determined by tu q in the trigonometric 
functions. For 5 = one finds exactly jg = 2rcos(g/2) which gives the uppermost curve in the 
figure when folded into the first Brillouin zone. 

Writing 7 max = vr — A, the maxima lead to a long-time behaviour of the correlations of the form 

cos(2n7 max ) = cos(2n(7r — A)) = cos(2nA) (29) 

and thus to oscillations with frequency 2A. The effect is the same as in Section ITTT1 for the case r 
near 7r/2 and may be called an Umklapp effect in the frequency space due to the discrete times. 
The double-peak structure in the dispersion can also lead to beats in S(t), if the two frequencies at 
the maximum and the minimum are close to each other. Such phenomena have also been observed 




Finally, we show in Fig. fTS] how the asymptotic value of S depends on the parameters r and 
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FIG. 11: Dispersion 7 g for r = 0.77T and several values of the dimerization. 

5. The calculation uses the expressions (|41[) in the Appendix for the elements of the correlation 
matrix which is then diagonalized. In all cases, S(oo) is proportional to L. In the fully dimerized 
case, it is a periodic function of r with period it and maxima at odd multiples of ir/2. The latter 
times still play a special role for partial dimerization because the maximal value of oo q is always 
1, but S(oo) has only steps there. In general, it decreases as one lowers 6, which is reasonable, as 
then the changes during the quench become smaller. For 5 = 0, there is no change at all and S 
has to be constant in time and non-extensive. One should note that the asymptotic values give 
no information on how fast they are reached. An extreme case is 5 = 1 and r near kn. Then the 
approach is very slow and if r equals kn exactly, S remains at its initial value apart from periodic 
excursions as noted in Section [1111 
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FIG. 12: Asymptotic value of the entanglement entropy for L = 40 as a function of the quench period for 
different dimcrizations. 
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V. SUMMARY AND CONCLUSION 



We have studied an integrable quantum chain under the influence of a periodic quench. Because 
one switches back and forth between the two equivalent broken-symmetry arrangements, the time 
evolution is on average that of a translationally invariant system. The entanglement was studied 
via the entropy S for a segment of L sites. This contains more information than the Q measure 
related to reduced single-site density matrices and used in [7f. We found an overall behaviour 
similar to that in a single critical quench, namely an increase of S(t) to a final value which is 
proportional to the size of the subsystem. On a finer scale, the increase is seen to result from the 
time evolution within the quench periods which in the simplest case gives a step structure, but 
becomes partly oscillatory for longer periods and washed out if the dimerization is weak. On the 
scale of many periods, there are in general oscillations related to the spectrum of the evolution 
operator. They make the step heights slightly non-uniform and can lead to a striking wave-like 
approach of S to its asymptotic value. A particularly interesting special case was found for 5 = 1 
and r = 7r/2, where the electrons move like massless relativistic particles. This leads to very simple 
formulae and gives an elementary realization of the picture developed by Calabrese and Cardy for 
the entanglement process. 

We have stressed the close connections to the equilibrium statistics of two-dimensional lattice 
models. These were also visible in the studies of the kicked Ising chain [f], Q, |J but, somewhat 
surprisingly, not pointed out explicitly. The equilibrium models provide a useful reference frame. 
Thus the concept of the Hamiltonian limit, developed for very anisotropic equilibrium systems, 
could directly be taken over. On the other hand, the unitary time evolution also introduces some 
differences. The correlations are typically oscillating and the discrete time intervals between the 
periods give rise to lattice effects in the frequency domain. The quench could be classified as critical 
for all values of the period. Any modification in the relative length of the half-periods or in the 
two dimerizations would lead to gaps in the energies u q and thus to a non-critical quench. The 
entanglement in this case would still show similar overall behaviour due to the quasiparticles in H 
which propagate with the maximum velocity. However, the various band edges in the single-particle 
dispersions will make the picture more complex in detail. Similarly, a finite quench velocity would 



complicate matters since in our case one is permanently crossing the quantum critical 
The entanglement in such a crossing has so far only been studied for single quenches 



aoint 5 = 0. 
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S is not extensive but saturates for large L at a value depending on the transition time. It might 
be interesting to extend the considerations to the periodic case. In the framework of equilibrium 



models, this would correspond to a lattice with slowly varying periodic layering. 

Note added: A similar situation as studied here was investigated for Heisenberg chains in [361 ] . 
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Appendix: Calculation of the correlation matrix 

In this Appendix we summarize the method of obtaining the matrix C(f) from which the 
entanglement entropy is calculated for arbitrary dimerizations and quench periods. 

In order to calculate the time evolution of the fermionic operators a,j and bj, it is convenient to 
switch to the Fourier-transformed operators defined in ([2|). Then the a q and b q evolve for each q 
independently as 




COS LJqt ie l P q sinu)qt 




v q (t) q (30) 



ie t!Pq sin u q t cos uj q t I \ b„ j \ b n 

which is the analogue of (|15p . In Fourier space, the full evolution matrix always decomposes into 
such blocks. 

During the quench the system is evolved with two different Hamiltonians Hq, Hi with corre- 
sponding matrices vo q (r), vi q (r) and the time dependent operators after n full periods are given 
through 

u q (2nr) = [vo q {T)v lq (T)] n . (31) 

In order to obtain u q {2nr) one has to raise the matrix u q {2r) to the n-th power, thus one has 
to determine its eigenvalues and eigenvectors, as well. The eigenvalues can again be written as 
exp(±z7 g ) where 

cos 7 9 = cos 2 u> q T — cos((/?i — (po) sin 2 Ul q T (32) 



with 



cos 2 (g/2)-J 2 sin 2 (g/2) 

cos^i - ipo) = . 2 . (33) 

cos z ((//2) + d z sin (<7/2) 



To simplify the notation we have dropped the q indices of <po and ip\. Denoting the first and second 
component of the z-th eigenvector by A{ and Bi, respectively, the equations then yield 



l^l,2| 2 



1=F 



sm(<£>i — ipo) sin 2 u> q T 



and 



B 



1,2 



sin u) q T cos uj q T (e^ 1 + e^° ) 



A 



1,2 • 



(34) 



(35) 



± sin7q — sin((^i — <^o) sin 2 w 9 r 
With the help of (|32p -(|35 p one can construct the matrix u q (2nr) = u q (2r) whose elements read 



un = cos n7„ — i sm((^i — <po j sm uj„t 

smj q 



sm 



U\2 



(36) 



smj q 



+ e llf0 ) Sin UJ q T COS LO q T 



while the other two elements u 2 \ = — u\ 2 and u 22 = u*i follow from unitarity. 

Finally, one has to calculate the correlation matrix elements Cij(2nr) as defined in (I12p which 
requires knowledge of the initial correlations. These are given as 



(aU q ) = (b%) = 1/2 , (ap q ) = (b\a q y = e^/2 



(37) 



Going over to Fourier transforms and using the form of the evolution matrix u q {2nr) one obtains 
Cij(2nr) = \{8ij + Tij) with the block matrix 



/ 



n n_i ••• n^iv 

Hi n : 



n, 




(38) 



\njv_i n j 

where m = I — j and the matrix elements of the 2x2 matrices H m read in the thermodynamic 
limit 



u 



9n 



v da 

^-Re(u 22 u 12 e l ^)e- l<im , 

-7T T 

^ (u 2 22 e^° - u 2 2l e-^°) e~ iqm 



(39) 



2vr 



x 21 c 



Up to now we have considered only the special time instances t n = 2nr. To extend the calcu- 
lation to arbitrary times t n < t < t n+ \ one has to calculate the matrix elements of 



u q (t) = < 



Vlq(t - t n ) U g (t n ), if t<t n + T 

^v£ q (t n+ i - t) u q (t n+1 ), iit>t n + T 



(40) 



and substitute them in (|39p to obtain the matrix T for intermediate times. 

It is also interesting to investigate the asymptotic (n — * oo) behaviour of the matrix elements. 
Then the integrals in (I39p have rapidly oscillating arguments, and one can replace cos 717,2 anc ^ 
sin 717,2 zero as weu as c °s 2 n 7q an d sin 2 777,2 by 1/2, as these factors are multiplied by smooth 
functions of q and integrated over. After a lengthy calculation one has the form 

n dq i sin oj q T cos u q T 



f m (oo)= — -pr- — 2~g— 2 5 tan - sin , 

7_,r 27r o l tan z § + cos z uj„t 2 

/•* J 2 ( 41 ) 

g m {oo) = 7Tl2l — — 2 (cosgm + tan-smgm). 

j.^ 2ir o l tan z | + cos" 1 cj g r 2 

For 5 = 1 and r — ► all expressions simplify considerably. Then f m and g m become Bessel 
functions and their asymptotic values are / m (oo) = and g m (oo) = (5 m fl +<5 m ,i)/2. This gives the 



correlation matrix found in 



31| for a single quench. 
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